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MASSACHUSETTS 


ABSTRACT 


A  multi-dimensional  sequence  is  not,  in  general,  uniquely  defined  in 
terms  of  only  the  phase  or  magnitude  of  its  Fourier  transform.  However,  in 
this  paper  some  conditions  are  developed  under  which  a  multi -dimensional 
sequence  is  uniquely  defined  by  its  phase.  A  similar  set  of  conditions  are 
then  developed  for  the  unique  specification  of  a  multi-dimensional  sequence 
in  terms  of  its  Fourier  transform  magnitude.  In  both  cases,  it  is  initially 
assumed  that  either  the  phase  or  magnitude  is  known  for  all  frequencies.  The 
results  are  then  generalized  to  the  case  in  which  the  phase  or  magnitude  is 
known  only  for  a  finite  set  of  frequency  values. 
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I .  INTRODUCTION 


Under  a  variety  of  conditions,  a  multi -dimensional  (m-D)  sequence  may  be 
reconstructed  from  partial  information  about  its  Fourier  transform.  For  example, 
if  an  m-0  sequence  x(nj, . . .  ,0^^)  is  zero  whenever  any  of  the  indices  for 
k=l,...,m  are  negative,  then  x(nj,. . . ,0^^)  can  be  exactly  recovered  from  the  real 
part  or,  except  for  x(0,...,0),  from  the  imaginary  part  of  Its  Fourier 

transform  [1].  If,  on  the  other  hand,  x(nj . n^^)  is  minimum  phase  [2],  then 

it  can  be  recovered  from  the  magnitude  c. ,  to  within  a  scale  factor,  from  the 
phase  of  its  Fourier  transform. 

The  reconstruction  of  an  m-D  sequence  from  such  partial  information  is 
important  and  useful  in  many  applications.  For  example,  in  some  cases  of 
optical  image  processing  or  in  the  measurement  of  diffraction  patterns,  only 
spectral  magnitude  information  can  be  recorded  or  is  available.  Therefore,  it 
is  of  interest  to  recover  a  signal  from  only  spectral  magnitude  information  in 
these  cases.  In  other  situations,  either  the  spectral  magnitude  or  phase  of  a 
signal  may  be  severely  distorted  so  that  the  restoration  of  the  signal  must  rely 
on  the  undistorted  component.  For  example,  in  the  class  of  problems  referred  to 
as  blind  deconvolution  [3|,  a  signal  is  to  be  recovered  from  an  observation  which 
is  the  convolution  of  a  desired  signal  with  some  unknown  signal.  Since  little 
is  usually  known  about  either  the  desired  signal  or  the  distorting  signal, 
deconvolution  of  the  two  signals  is  generally  a  very  difficult  problem.  However, 
in  the  special  case  in  which  the  distorting  signal  is  known  to  have  a  Fourier 
transform  which  is  purely  real,  the  tangent  of  the  phase  of  the  observed  signal 
is  identical  to  the  tangent  of  the  phase  of  the  original  signal.  Such  a  situation 
occurs,  at  least  approximately,  in  long-term  exposure  to  atmospheric  turbulence 
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or  when  images  are  blurred  by  severely  defocused  lenses  with  circular  aperture 
stops  [4).  In  this  case,  it  is  of  interest  to  consider  signal  reconstruction 
from  phase  information  alone. 

This  paper  is  concerned  with  the  development  of  some  conditions  under  which 
an  m-D  sequence  is  uniquely  defined  by  its  Fourier  transform  phase  or  magnitude. 
In  general,  of  course,  phase  or  magnitude  information  alone  is  not  sufficient 
to  uniquely  specify  an  m-U  sequence.  For  example,  any  m-D  sequence  may  be 
convolved  with  a  zero  phase  sequence  to  produce  another  m-0  sequence  with  the 
same  phase.  Similiarly,  any  m-D  sequence  may  be  convolved  with  an  all-pass 
sequence  to  produce  another  m-D  sequence  with  the  same  magnitude.  Therefore, 
without  any  additional  information  or  constraints,  the  Fourier  transform  phase 
or  magnitude  may,  at  best,  uniquely  specify  an  m-D  sequence  to  within  a  zero 
phase  or  all-pass  convolutional  factor,  respectively.  Nevertheless,  with  a 
few  basic  results  from  the  theory  of  polynomials  in  several  variables,  some 
useful  conditions  may  be  derived  under  which  an  m-D  sequence  is  uniquely  defined 
by  the  phase  or  magnitude  of  its  Fourier  transform.  These  conditions,  which  are 
distinctly  different  from  the  minimum  or  maximum  phase  constraints,  imply  that 
most  m-D  sequences  with  finite  support  are  recoverable  from  either  their  phase 
or  magnitude. 

This  paper  is  organized  as  follows.  In  Section  II,  the  necessary  results 
from  the  algebra  of  polynomials  in  more  than  one  variable  are  briefly  reviewed. 
Some  notation  and  terminology  related  to  multi-dimensional  signals  is  then 
presented  in  Section  III.  In  Section  IV,  conditions  are  developed  under  which 
an  m-D  sequence  is  uniquely  defined  in  terms  of  the  phase  or  magnitude  of  its 
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Fourier  transform.  Although  these  conditions  are  initially  derived  under  the 
assumption  that  either  the  phase  or  magnitude  is  known  for  all  frequencies, 
they  are  then  extended  to  the  case  in  which  the  phase  or  magnitude 
is  known  only  over  a  finite  set  (lattice)  of  points. 


II.  POLYNOMIALS  IN  MORE  THAN  ONE  VARIABLE 


In  this  section,  some  notation  and  terminology  related  to  the  algebra  of 
polynomials  in  more  than  one  variable  is  reviewed.  In  addition,  two  theorems 
are  presented  which  are  of  considerable  importance  in  many  multi-dimensional 
signal  processing  applications  and  will  be  referred  to  frequently  in  this  paper. 
Proofs  of  these  theorems  as  well  as  a  detailed  treatment  of  many  topics  not 
presented  in  this  section  may  be  found  in  [5]. 

A  monomial  is  a  function  of  n  variables  of  the  form: 


(1) 


where  kj,k2,. . . ,k^  are  non-negative  integers  and  c  is  an  arbitrary  number  which 
is  referred  to  as  the  coefficient  of  the  monomial.  The  degree  of  the  monomial 
is  defined  as 


d(f)  =  kj  +  k2  +  ...  +  k^  (2) 

A  polynomial  is  then  simply  a  sum  of  a  finite  number  of  monomials.  In  this 
paper,  a  polynomial  will  be  referred  to  as  being  non-trivial  if  it  consists  of 
a  sum  of  two  or  more  monomials.  A  trivial  polynomial  is  therefore  either  a 
constant  or  a  monomial  of  non-zero  degree. 

The  coefficients  of  the  monomials  which  define  a  polynomial  are  called  the 
coefficients  of  the  polynomial  and  the  degree  of  the  polynomial  is  the  degree 
of  the  monomial  with  the  highest  degree.  Any  polynomial  p  in  n  variables  of 
degree  N  may  therefore  be  written  in  the  form: 
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(3) 


kj+. .  .+k^£N 


It  is  often  useful  to  consider  p  In  (3)  as  a  polynomial  in  one  variable, 
say  with  coefficients  which  are  polynomials  in  the  remaining  (n-1)  variables. 
For  example,  p  in  (3)  may  be  written  as 

N 

k=0 

where  (})o(k)  for  k=0,l,.,.,M  are  polynomials  in  the  (n-1)  variables  z  for  n^£. 
x-  n  ' 

In  this  form,  the  largest  value  of  k  for  which  <i>j^(k)  is  non-zero  is  referred 
to  as  the  degree  of  p  with  respect  to  the  variable  and  will  be  denoted  by 
clj(p). 

If  all  the  coefficients  of  a  polynomial  p  belong  to  a  particular  number 
field,  ^  then  p  is  called  a  polynomial  over  the  field  3^.  The  set  of  all 
polynomials  in  n  variables  over  a  field  ^ form  a  ring  which  is  denoted  by 

^(Zp..,,z^).  If  two  polynomials  p^  and  p^  in  ^(Zj . z^)  are  such  that 

Pj=cpj  where  ce,^ and  is  non-zero,  then  p^  and  called  associated 

polynomials.  A  polynomial  p  e  jr(Zj,. . .  ,z^)  with  d(p)>0  is  called  a  reducible 
polynomial  over  the  field  ^ if  there  exists  polynomials  Pj,P2  •  •  iZ^) 

such  that  p=PjP2  with  d(pj)>0  and  d(p2)>0.  If  no  such  decomposition  is 
possible,  then  p  is  called  an  irreducible  polynomial.  It  is  of  interest  to 
note  that,  as  a  consequence  of  the  Fundamental  Theorem  of  Algebra,  the  only 
polynomials  in  one  variable  over  the  field  of  complex  numbers  which  are 
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irreducible  are  polynomials  of  first  degree.  This  result,  however,  is  not  true 
in  general.  Over  the  field  of  real  numbers,  for  example,  the  polynomial 
p(z)=z  +1  is  irreducible.  Even  over  the  field  of  complex  numbers,  the  polynomial 
p(Zj,Z2)=z^+Z2  is  irreducible.  This  second  example  illustrates,  in  particular, 
the  fact  that  the  Fundamental  Theorem  of  Algebra  does  not  hold  for  polynomials 
in  more  than  one  variable. 

The  first  theorem  of  interest  in  this  paper  asserts  that  any  polynomial  of 
non-zero  degree  can  be  uniquely  decomposed,  to  within  factors  of  zero  degree, 
into  a  product  of  irreducible  polynomials.  More  specifically. 


Theorem  1:  Any  polynomial  pe^(z^,...  ,z^)  having  non-zero  degree 
can  be  expressed  as  a  product  of  factors  irreducible  in 
Furthermore,  if  p  has  two  different  factorizations: 


P  =  =  91  9, 


then  m=n  and  the  factors  f^-  and  g^  can  be  ordered  in  such  a  way 
that  the  factors  are  associated. 


It  is  well-known  [5]  that  a  polynomial  p(z)  in  one  variable  of  degree  N 
is  uniquely  defined  in  terms  of  its  values  over  a  set 

distinct  points.  This  result  has  a  natural  extension  to  polynomials  in  more 


than  one  variable  which  asserts  the  uniqueness  of  a  polynomial  p  e^{Zj, . . .  ,z^) 
in  tenns  of  its  values  over  an  n-dimensional  lattice  of  points.  An  n-dimensional 
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lattice  is  an  extension  of  a  single  set  of  points.  A,  to  an  n-fold  Cartesian 
product  of  n  sets  of  points,  for  k=l,...,n.  Specifically,  let  A|^  =  |a|^ 
be  a  set  of  Nj^  distinct  points  in  the  field  ^  for  k=l,...,n.  Then  the 
n-dimensional  lattice  A j,. .. ,A^)  is  defined  as 

n 

^(Ap....A^)  =  n  =  A^  X  A2  X  ...  X  A^  (6) 

k=l 

The  result  of  interest  may  now  be  stated  as  follows: 

Theorem  2:  Suppose  p^,P2  c^lzj,...  ,z^)  with  d|^(pp  <  Nj^  and 
dk(p2)  <  for  k=l,...,n.  Let  A,^  be  a  set  of  N,^  distinct  points 
in  the  field  If  p^  and  P2  are  equal  over  the  set  (lattice)  of 
points  i?(Aj,...,A^),  then  Pi=P2- 

This  theorem  will  be  used  in  deriving  some  conditions  for  the  uniqueness  of  a 
multi-dimensional  sequence  in  terms  of  a  finite  set  (lattice)  of  values  of 
either  the  phase  or  magnitude  of  its  Fourier  transform. 
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III.  FRAMEWORK 


This  paper  is  concerned  with  the  uniqueness  of  a  multi -dimensional  sequence 
with  real  coefficients  in  terms  of  either  the  phase  or  magnitude  of  its  Fourier 
transform.  Although  the  results  which  are  presented  apply  to  sequences  of 
arbitrary  dimension,  all  discussions  are  phrased  in  terms  of  two-dimensional 
sequences,  x(nj,n2),  in  order  to  simplify  notation.  In  this  section,  some 
notation  and  terminology  related  to  2-D  sequences  is  presented  and  the  general 
framework  of  the  uniqueness  problem  is  established. 

The  z-transform  of  a  two-dimensional  sequence,  x(nj,n2)  is  defined  by 


n.  n2 

x(nj,n2)  Zj  Z2 


(7) 


"1  "2 

and  the  Fourier  transform,  denoted  by  X(a)j,aj2)»  equal  to  X(Zj^,Z2)  on  the  unit 
bi-disc  jZjl  =  jZ2|  =  1,  i.e.. 


X(a)j,0i)2)  ~  X(Zj^,Z2) 


Zj=exp(ja,j) 

Z2=exp(ju)2) 


(8) 


Written  in  polar  form,  X(aj2,w2)  is  represented  in  terms  of  its  magnitude  and  phase 
as 

X((i)j^  ,^2)  ~  I  »(i)2 )  1  j *1*^2 


where  it  is  assumed  that  the  phase  is  defined  by  its  principal  value.  Therefore, 
in  terms  of  (9),  this  paper  considers  the  uniqueness  of  x(nj,n2)  in  terms  of 
(|)j^((Dj,(D2)  or  |X(a)j ,0)2)1  •  Although  initially  it  is  assumed  that  either  the  phase 
or  magnitude  is  known  for  all  values  of  o)j  and  0)2,  the  results  are  then  extended 
to  the  case  in  which  the  phase  or  magnitude  is  known  only  over  a  finite  set  of  points. 
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Most  of  the  sequences  which  are  considered  in  this  paper  have  finite 
support,  i .e.,x(n^,n2)  is  non-zero  only  for  finitely  many  values  of  its  arguments 
n^^  and  n^.  A  sequence  with  finite  support  is  said  to  be  of  extent  NjXN2  if 
x(nj,n2)=0  outside  a  rectangular  region  of  the  form  x  CLq»Lq+N2-1] 

for  some  and  L^.  For  the  special  case  in  which  and  may  be  taken  to  be 
equal  to  zero,  x(n^,n2)  is  said  to  have  first  quadrant  support  and  the  region  of 
support  is  denoted  by  *  i-e., 

=  [O.Nj-1]  x  [0,N2-1]  (10) 

Throughout  this  paper,  any  sequence  with  finite  support  may  be  assumed,  without 
any  loss  in  generality,  to  have  first  quadrant  support.  In  the  general  case,  a 
sequence  may  simply  be  shifted  in  order  to  satisfy  this  assumption.  Consequently, 
the  set  F(n^,n2)  will  be  used  to  denote  the  collection  of  all  real  2-D  sequences 
of  finite  extent  with  first  quadrant  support.  The  notation  XEF(nj,n2)  will 
therefore  mean  that  the  sequence  x(nj,n2)  has  support  ^(^■^,^2)  for  some  Nj  and 
N2.  Since  the  z-transform  of  a  sequence  XEF(nj,n2)  is  a  polynomial  in  two 
variables,  z^  and  Z2,  X(zj,Z2)  is  an  element  of  ^(Zj,Z2)  over  the  field  of 
real  numbers: 

xeF(nj,n2)  Xe  ^^(Z'^tZ2)  (11) 

Now,  suppose  that  xcFCn^.ng)  and  has  an  irreducible  z-transform  of  degree 
Nj  in  Zj  and  N2  in  Z2  and  consider  the  sequence  xeFCn^.ng)  defined  by 
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x(n^,n2)  =  x(  N^-n^,  N^-n^)  (12) 

The  z-transform  of  x(nj,n2)  is  also  irreducible  and  is  given  by 

-1  -1 

X(Zj,Z2)  =  Zj  Z2^  X(Zj  ,Z2^)  (13) 

A  special  class  of  sequences  which  will  play  an  important  role  in  the  following 
sections  consists  of  those  sequences  for  which 


x(nj,n2)  =  ±x(nj,n2) 

(14) 

X(Zj,Z2)  =  ±X(Zj,Z2) 

(15) 

Since  '^(0^,02)  corresponds  to  a  180®  rotation  of  x(nj,n2)»  sequences  which 
satisfy  (14)  are,  except  possibly  for  a  minus  sign,  invariant  under  180° 
rotations.  Therefore,  these  sequences  will  be  said  to  be  symmetric  or  to  have 
symmetric  z- transforms.  It  should  be  noted  that  a  1-D  sequence  which  has  a 
symmetric  z-transform  has  all  of  its  zeros  (excluding  those  at  z=0  or  z'^=0) 
on  the  unit  circle  or  in  reciprocal  pairs.  Therefore,  (14)  represents  an 
extension  of  this  property  to  2-D  sequences.  It  may  also  be  noted  that  except 
for  a  linear  phase  term  of  the  form  exp[Jsj(N2a)j+N2U)2)] ,  the  Fourier  transform 
of  a  symmetric  sequence  is  either  purely  real  or  purely  imaginary. 

An  important  property  of  a  sequence  xeF(nj,n2)  is  that  its  z-transform 
need  only  be  known  over  a  finite  set  of  points  in  order  to  uniquely  specify 
the  sequence.  Although  these  points  can  not  be  chosen  arbitrarily.  Theorem  2 
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in  Section  II  provides  one  set  of  points  which  is  sufficient  for  this 
specification.  More  specifically,  the  following  lemma  is  a  direct  consequence 
of  Theorem  2: 


Lemma:  Suppose  x,yeF(nj,n2)  with  support  ^(Nj,N2).  Let 
^k"{^k  j}  ®  *^k  complex  numbers  for  k=l,2 

with  >  Nj  and  M2  >  N2.  If 


X(Z,,Zp)  =  Y(z,,Z2) 

'  ^  ^(A,.A2)  *  ^  ^’(Aj.Aj) 

then  x(nj,n2)  =y(nj,n2)  for  all  n^  and  02* 


Of  particular  interest  in  this  paper,  however,  is  the  case  in  which  the 
elements  of  the  sets  A|^  have  unit  magnitude.  Specifically,  consider  the  sets 


(17a) 


(17b) 


where  the  elements  of  for  k=l,2  are  assumed  to  be  distinct.  Then 


X(z.,zJ  =  X(w.,a)J(  (18) 

^  ^  <2'(Ai,A2)  ^ 

represents  the  Fourier  transform  of  x(nj,n2)  evaluated  over  the  lattice 
SC  the  ajj,w2-plane.  If,  in  addition,  the  numbers  ^  are  equally 
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spaced  between  0  and  2tt,  j^=2Tril/Mj^,  then  (18)  represents  the  MjXMg-point 
Discrete  Fourier  Jransform  (DFT)  of  x(npn2).  In  this  case,  the  MjXM2-point  DFT 
will  be  denoted  by  X(kj,k2)|^^^Hi|^: 


X(k2,k2)|^ 


l’^2 


X(Zj,Z2) 


Zj=exp(j2Trkj/Mj) 

Z2=exp(j2TTk2/M2) 


(19) 


which,  when  expressed  in  terms  of  its  magnitude  and  phase,  will  be  written  as 


X(ki,k2)f^^^M2  ^ 

In  the  following  sections,  when  the  length  of  the  DFT  is  assumed  to  be  known 
or  when  it  is  explicitly  stated,  the  subscripts  in  (20)  will  occasionally  be 
dropped. 
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IV.  THE  QUESTION  OF  UNIQUENESS 

In  this  section,  the  results  from  Section  II  are  applied  to  the  question 
of  the  uniqueness  of  a  2-D  sequence  in  terms  of  the  phase  or  magnitude  of  its 
Fourier  transform.  In  Section  4.1,  conditions  are  given  under  which  a  2-D 
sequence  with  finite  support  is  uniquely  specified  in  terms  of  the  phase  of  its 
Fourier  transform.  A  similiar  set  of  conditions  are  given  in  Section  4.2  for 
the  uniqueness  of  a  2-D  sequence  with  finite  support  in  terms  of  its  Fourier 
transform  magnitude.  In  Section  4.3,  the  results  in  Sections  4.1  and  4.2  are 
used  to  generate  a  set  of  conditions  under  which  a  2-D  sequence  with  finite 
support  in  uniquely  specified  in  terms  of  either  its  Fourier  transform  phase 
or  magnitude.  In  addition,  the  generalization  of  these  results  to  2-D 
sequences  whose  convolutional  inverses  have  finite  support  is  described. 

4.1  Uniqueness  in  terms  of  phase 

It  has  recently  been  shown  [6]  that  a  1-D  finite  length  sequence  x(n)  is 
uniquely  specified  to  within  a  scale  factor  by  the  phase  or  the  tangent  of  the 
phase  of  its  Fourier  transform  if  X(z)  has  no  zeros  on  the  unit  circle  or  in 
reciprocal  pairs,  i.e.,  if  X(z)  contains  no  symmetric  factors.  This  result  may 
be  directly  extended  to  the  case  of  multi-dimensional  sequences.  This  extension, 
in  terms  of  2-D  sequences  is  as  follows: 

Theorem  3:  Let  x,y e  F{nj,n2).  If  X(Zj,Z2)  and  Y(z^,Z2)  have  no 
symmetric  factors  and 

*^x^“r“2^  =  <l>y(wj,w2)  (21a) 
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for  all  and  Wg,  then  y{nj,n2)  =  ex(nj,n2)  for  some  positive 
number  3.  If.  on  the  other  hand, 

tan  =  tan  (j)y(a)j,ti)2)  (21b) 

for  all  (jjj  and  then  y(nj,n2)  =  3x(nj,n2)  for  some  real  number  3. 

It  should  be  noted  that  the  symmetric  factors  which  are  excluded  from  X(Zj,Z2) 
and  Y(Zj,Z2)  in  this  theorem  need  not  be  irreducible.  For  example,  if  A(Zj,Z2) 

'\f 

is  a  polynomial  and  X(Zj,Z2)=P(2j, 22)0(2^.22)  where  P(Zj,Z2)=A(Zj,Z2)A(2j,Z2) , 
then  x(nj,n2)  does  not  satisfy  the  constraints  of  the  theorem  since  P(Zj,Z2)  is 
a  (reducible)  symmetric  factor  of  X(zj,Z2).  In  effect,  the  exclusion  of  symmetric 
factors  from  a  sequence  xeF(nj,n2)  is  equivalent  to  the  constraint  that  if 
A(Zj,Z2)  is  an  irreducible  factor  of  X(Zj^,Z2)  then  A(z^,Z2)  is  not  a  factor  of 
X(  z  j ,Z2 )  • 

An  outline  of  a  proof  of  Theorem  3  is  as  follows.  Let  x,yeF(nj,n2)  and 
let  N  be  a  positive  integer  which  is  sufficiently  large  so  that  x(nj,n2)  and 
y(nj^,n2)  are  both  zero  outside  the  domain ^(N,N) .  Consider  the  sequence 

h(nj,n2)  =  x(nj,n2)  *  y(-np-n2)  (22) 

which  has  a  z-transform  given  by 


H(z,,z,)  =  Y(z7^,Z9^) 


(23) 


By  noting  that  the  phase  of  hCn^^.n^)  is  given  by 


(24) 


it  follows  from  {21a)  that  4)^(ojj,W2)=0  or  from  (21b)  and  the  trigonometric 
identity 

tan  a  +  tan  B 

tan(a+B)  =  - 

1  +  tan  a  tan  6 

that  tan(t.f,(ujj,a)2)  =0.  In  either  case,  the  Fourier  transform  of  h(nj,n2)  is 
purely  real  which  implies  that 


h(nj,n2)  =  hC-n^.-ng) 


(26) 


Therefore,  from  (23)  and  (26), 


X(Zj,Z2)  Uzl^alh  =  X(z’^Z2^)  Y(Zj,Z2)  (27) 

Multiplying  both  sides  of  (27)  by  (ZiZ2)'^  results  in  the  following  polynomial 
equation: 


^  k,  kp  'v.  \  ^1 

X(Zj,Z2)  Y(z^,Z2)  z^^  z/  =  XlZj,Z2)  Y(Zj.Z2)  z^  Z2 


(28) 


where  k2,k2,tj,  and  I2  non-negative  integers.  Now  consider  an  arbitrary 
non-trivial  irreducible  factor  X|^(Zj,Z2)  of  X(Zj^,Z2).  From  Theorem  1  in 
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Section  2,  it  follows  that  X,^( 2^,22)  must  be  associated  either  with  a  factor  of 
X(2^,22)  or  with  a  factor  of  Y(2j,22}.  However,  if  associated  with 

a  factor  of  Xj^(2j,22)  then 

X^(2i,22)  =  aXj,(2j,22)  (29) 

for  some  1.  If  a=k  then  (29)  implies  that 

=  Qt^  X,^(2j.22)  (30) 

Therefore,  a  =  il  and  Xi((ziiZ2)  is  s^etric.  If,  on  the  other  hand,  k?**,  then 

Xk(2j,22)  X^(2j,22)  =  aXjj(2j,22)  Xj^(2i,22)  (31) 

In  other  words,  both  Xj^(2j,22)  and  Xj^(2j,22)  are  factors  of  X{2j,22).  Both 
cases,  however,  are  excluded  by  the  theorem  hypothesis.  Consequently,  each 
non-trivial  irreducible  factor  of  Xltj.tg)  must  be  associated  with  a  factor 
of  Y(2j,22).  Therefore,  Y(2j,22)  is  of  the  form: 

tm,  m,  ■)  ,  , 

2j  22^  P(2i,Z2)J  X(2i,22)  (32) 

where  m^  and  m2  are  integers  and  P(2j,22)  is  a  polynomial.  However,  as  in  the 

steps  leading  to  (26),  (21)  requires  that  the  term  in  brackets  correspond  to  a 
sequence  with  a  real  Fourier  transform.  This,  in  turn,  implies  that  P(2j,22) 
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is  symmetric  which,  since  Y(Zj,Z2)  contains  no  symmetric  factors,  requires  that 
P(Zj,Z2)  be  a  constant; 


m.  m^ 

Y(2j,Z2)  *  6Zj  Z2  X(Zj,Z2) 


Again  using  (21)  it  follows  that  mj=m2=0  and  the  desired  result  follows  by 
noting  that  6  must  be  positive  in  the  case  of  hypothesis  (21a). 

It  may  be  apparent  at  this  point  that  it  is  not  necessary  to  assume  that 
(21)  holds  for  all  and  0)2  in  order  to  prove  the  theorem  if  the  regions  of 
support  of  x(n^,n2)  and  y(n^,n2)  are  known.  More  specifically,  suppose  that 
x(n2,n2)  and  y(nj,n2)  are  known  to  be  zero  outside  the  domain  ^(N^,N2). 

If  Mj>2(Nj-l)  and  M2>2(N2-1)  and  if  (21)  is  replaced  with  the  constraint 
that  the  phase  or  tangent  of  the  phase  of  the  M2xM2-point  DFT’s  of  x(nj,n2) 
and  y(nj,n2)  are  equal,  then  Theoren  3  still  remains  valid.  This  follows  first 
from  the  observation  that  h(nj,n2)  in  (22)  is  zero  outside  the  region 
[(1-Nj) ,(Nj-1)]x[(1-N2) ,(N2-1)]  .  Therefore,  the  MjXM2-point  DFT  of  h(nj,n2) 
equals  the  product  of  the  MjXM2-point  DFT's  of  x(nj,n2)  and  y(nj,n2).  Thus, 


Mj,M2 


Mj,M2 


Mj,M2 


and  (27)  follows  as  in  the  proof  of  Theorem  3.  Proceeding  as  in  Theorem  3,  it 
follows  that  Y(Zj,Z2)  and  X(zj,Z2)  are  related  by  (32).  However,  the  constraints 
on  the  phase  of  the  MjXM2-point  UFT's  of  x(nj,n2)  and  y(n^,n2)  imply  that  the 
term  in  brackets  in  (32)  corresponds  to  a  sequence  with  a  real  MjxM2>point  DFT. 
Therefore,  in  a  style  similiar  to  that  in  Theorem  3,  it  follows  that  the  term 
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in  brackets  in  a  constant.  Therefore  a  corollary  to  Theorem  3  is  as  follows: 


1 


Corollary  3.1:  Let  x,yeF(nj,n2)  with  support  i^(Nj,N2).  If 
X(Zj,Z2)  and  Y(zj,Z2)  have  no  symmetric  factors  and 


<l)y(ki,k2)  =  <t>x(ki.k2^ 


Mj,M2 


Hj,M2 


(35) 


with  Mj>2(N^-l)  and  M2>2(N2-1).  then  y(nj,n2)=  B  x(nj,n2)  for  some 
positive  number  6.  If,  on  the  other  hand, 


tan  (t)„{k,  .k,)  =  tan  4>^(k,,kp) 

^  ^  2  Mj,M2  ^  ^  ^  Mj,M2 

then  y(nj,n2)* B x(nj,n2)  for  some  real  number  B. 


(36) 


The  importance  of  this  corollary  lies  in  the  fact  that  it  allows  for  the 
development  of  practical  algorithms  for  reconstructing  a  sequence  from  the  phase 
of  its  DFT. 

4.2  Uniqueness  in  terms  of  magnitude 

In  Section  4.1,  the  uniqueness  of  a  multi-dimensional  sequence  in  terms  of 
the  phase  of  its  Fourier  transform  was  considered.  This  section  addresses  the 
dual  problem  related  to  the  uniqueness  of  a  multi -dimensional  sequence  in  terms 
of  its  Fourier  transform  magnitude.  Perhaps  the  first  treatment  of  this  question 
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of  uniqueness  may  be  found  in  (7l  where  the  uniqueness  of  a  2-D  sequence 
xeF(nj,n2)  is  stated  to  be  related  to  the  irreducibility  of  its  z-transform. 

In  this  section,  a  slightly  more  general  result  is  derived  which  includes 
sequences  with  irreducible  z-transforms  as  a  special  case.  Even  more  importantly, 
however,  as  in  section  4.1,  the  uniqueness  of  a  sequence  xeF{nj,n2)  is  also 
considered  when  the  magnitude  of  its  Fourier  transform  is  known  only  over  a 
finite  set  (lattice)  of  points. 

Consider  a  sequence  xeF(nj,n2)  for  which  |X(a)j,(i)2)  1  is  known  for  all 
and  W2-  Since  the  inverse  Fourier  transform  of  |X(a)j^,t02) is  the  auto¬ 
correlation,  x(nj,n2): 

rj^(n2,n2)  -  x(n2,n2)  *  x(~n2,-n2)  (37) 

the  specification  of  |X(a)2,(jD2)  1  is  equivalent  to  the  knowledge  of 
its  z-transform; 

Rx(Zi.Z2^  =  X(Z2,Z2)  X(z'^,Z2^)  (38) 

For  any  xeF(n2,n2),  the  most  general  form  for  its  z-transform,  X(Z2,Z2) 
is  given  by 

m,  m, 

X(Z2,Z2)  =  Z2^  n  ^k^^l’^2^ 

k=l 

where  8  is  a  real  number,  mj  and  m2  are  non-negative  integers,  and  X|^(z2,22) 
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for  k=l,...,p  are  non-trivial  irreducible  polynomials.  Substituting  (39)  into 


(38)  gives 


Rx(Zi,Z2) 


P 

^  n  ^k^^l  *^2  ^ 

k=l 


(40) 


Now  suppose  that  the  polynomial 


P 

P(Zi.Z2)  =  n 
k=l 

N,  N^ 

is  of  degree  in  and  N2  in  Z2.  Multiplying  R^(Zj,Z2)  by  Zj  Z2  yields  a 
polynomial  in  and  Z2  which  has  degree  2Nj  in  z^  and  2N2  in  Z2: 

N,  N2  2  ^  ''' 

V  ^x^^l’^2^  =  P  n  X,^(Zj,Z2)  X,^(Zj,Z2)  (42) 

k=l 


It  is  apparent  that  the  polynomials  Qx(zt,Z2)  and  |X(u)pa)2)|  contain  exactly 
the  same  information  about  x(nj,n2)  since  one  may  be  uniquely  derived  from  the 
other.  Therefore,  the  ability  to  uniquely  recover  x(nj,n2)  from  |X(u)j,a)2)l  is 
equivalent  to  the  ability  to  uniquely  recover  X(Zj,Z2)  from  Q^(;^j,Z2).  With  this 
in  mind,  it  follows  that  x(nj,n2)  cannot  be  unambiguously  recovered  from  magnitude 
information  alone.  For  example,  the  sign  of  B  as  well  as  the  linear  phase  term 

nip 

Zj  Z2^  is  not  recoverable  from  Qx^2i’Z2).  Even  more  importantly,  however,  is 
the  observation  that,  without  additional  information,  it  is  not  possible  to 
determine  whether  X|j(Zj,Z2)  or  X|^(zj,Z2)  is  a  factor  of  X(zj,Z2).  This  ambiguity 
is  not  surprising,  however,  since  it  represents  a  2-D  extension  of  a  familiar 
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result  for  1-D  sequences  |8l.  Specifically,  for  any  finite  duration  sequence 
x(n),  another  sequence  y(n)  may  be  generated  which  has  the  same  Fourier  transform 
magnitude  as  x(n)  by  simply  reflecting  a  zero  of  X(z)  about  the  unit  circle. 

For  2-D  sequences,  represents  the  reflection  of  the  zero  contour  of 

about  the  unit  bi-disc  |Zj|  =  \Z2\  ~  1* 

It  will  be  useful  in  the  following  discussions  to  define  an  equivalence 
relation  on  the  set  FCn^.n^)  as  follows: 

y(n^,n2)  'v<  x(nj,n2)  if  y(n^ ,02)  =  +  x(  n^ ,  k2+ n2  )  (43) 

for  some  integers  k^  and  k2.  In  other  words,  the  equivalence  class  generated 
by  a  sequence  x(nj,n2)  is  defined  to  be  the  set  of  all  sequences  which  may  be 
derived  from  x(nj,n2)  by  a  linear  shift,  a  time-reversal,  or  by  a  change  in 
sign  of  the  sequence.  Note  that  all  of  the  sequences  within  a  given  equivalence 
class  have  the  same  Fourier  transform  magnitude.  Thus,  it  will  be  convenient  to 
refer  to  the  Fourier  transform  magnitude  of  the  sequences  within  an  equivalence 
class  as  the  Fourier  transform  magnitude  of  the  class. 

In  general,  there  will  be  more  than  one  equivalence  class  having  the  same 
Fourier  transform  magnitude.  More  specifically,  given  a  sequence  xeF(nj^,n2) 
there  may  exist  another  sequence  yeF(nj^,n2)  with  the  same  Fourier  transform 
magnitude  as  x(nj,n2)  but  which  is  not  in  the  same  equivalence  class  as  x(nj,n2). 
Therefore,  the  goal  of  this  section  is  to  determine  a  set  of  conditions  which 
guarantee  the  existence  of  only  one  equivalence  class  with  a  given  Fourier 
transform  magnitude.  The  first  question  to  be  addressed,  however,  concerns  the 


number  of  equivalence  classes  which  have  a  given  Fourier  transform  magnitude. 
Once  this  has  been  established,  conditions  which  guarantee  the  existence  of  only 
one  equivalence  class  may  easily  be  determined.  The  answer  to  this  first 
question  is  implied  by  the  following  theorem: 

Theorem  4:  Let  xeF(nj,n2)  with  a  z-transform  given  by 
m,  m,  P 

X(z^,Z2)  =  ez^^  Z2  n  Xj^(z^,Z2) 

k=l 

where  ^ire  non-trivial  irreducible  polynomials  for 

k=l,...,p.  If  y e  F(nj,n2)  and 

iX(a)^,a)2)|  =  |Y(a)j,a)2)| 

for  all  ojj  and  u)2»  then  Y(zj,Z2)  is  of  the  form: 

£  il 

kel  k^I 

where  I  is  a  subset  of  the  integers  in  the  interval  [l,pj. 

This  theorem  is  simply  a  statement  of  the  fact  that  the  only  way  to 
generate  a  new  sequence,  y(nj,n2),  which  has  the  same  Fourier  transform 
magnitude  as  x(n^,n2)  is  to  convolve  x(nj,n2)  with  an  all-pass  sequence, 
h(nj,n2),  i.e.,  a  sequence  with  |H(a)j,u)2) i  =  1  for  all  and  However, 
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any  all-pass  sequence  with  a  rational  z-transform  is  of  the  form: 


H(Zj,Z2)  =  +  n 

k=l 

Therefore,  given  a  sequence  xeF(nj,n2)  with  a  z-transform  of  the  form  (44), 
y(n^,n2)  = x(nj,n2)  * h(nj,n2)  has  finite  support  if  and  only  if  for  each  k, 
Hk(Zi,Z2)  =  Xj^(zi,Z2)  for  some  ile[l,p].  Consequently,  Y(Zj,Z2)  must  be  of  the 
form  given  by  (46). 

An  outline  of  a  proof  of  this  theorem  is  as  follows.  With  x,y e F(nj,n2) , 
let  N  be  an  integer  large  enough  so  that  x(nj,n2)  and  y(nj,n2)  are  zero  outside 
the  domain  ^(N,N).  From  (45)  it  follows  that 


X(Zj,Z2)  X(zJ^Z2M  =  Y(Zj,Z2)  Y(zJ\z2^) 


-1 


(48) 


Therefore,  let  the  z-transform  of  y(nj,n2)  be  given  by 

z  z  ^ 

Y(Zj,Z2)  =  aZj^  Z2^  JJ  Y|^(Zj,Z2) 

k=l 


(49) 


where  Y|^(Zj^,Z2)  are  non-trivial  irreducible  factors  for  k=l,...,q.  Substituting 
(44)  and  (49)  into  (48)  and  multiplying  by  (21^2^^  yields  the  following  equation: 

M  M  ^  %  L  L  ^ 

6^Zj^Z2^nX^(z^,Z2)X^(Zj,Z2)  =  a^Zj^Z2^  nY,^(Zj,Z2)Y,^(Zj,Z2)  (50) 

k=l  k*»l 

where  Mj,  ^^2*  Lj,  and  L2  are  positive  integers.  From  Theorem  1  in  Section  2,  it 


23 


follows  that 

p  ^  q  ^ 

n  Xk(2i.Z2)  X,^(z^,Z2)  =  ct^  n  Y|^{Zj,Z2)  Y,^(Zj,Z2)  (51) 

k=l  k=l 

Again  from  Theorem  1,  it  follows  that  the  factors  Y|^(Zj,Z2)  may  be  ordered  in 

'X/ 

such  a  way  that  Y|^(Zj,Z2)  is  associated  with  either  X|^(Zj,Z2)  or  X|^(Zj,Z2)  for 
each  k.  Therefore,  from  (49)  and  the  fact  that  (45)  implies  a  =  ±g,  the  desired 
result  (46)  follows. 

It  should  be  noted  that,  as  a  consequence  of  this  theorem,  all  sequences 

in  F(nj,n2)  with  a  given  magnitude,  |X{u)j,w2)l.  have  z-transforms  with  the  same 

number,  p,  of  non-trivial  irreducible  factors.  Furthermore,  except  for  a  scale 

factor  of  (-1)  and  linear  shifts,  the  only  way  to  generate  another  sequence 

y(n2>n2)  in  F(n2,n2)  with  the  same  magnitude  as  x(nj,n2)  is  to  replace  one  or 

more  non-trivial  factors  X|^(Z2,Z2)  of  X(z2,Z2)  with  X|^(Zj,Z2).  However,  if 

''' 

is  symmetric,  then  the  replacement  of  X|^(Z2,Z2)  with  X|^(Z2.Z2)  may 
only  change  X(Z2,Z2)  by  a  factor  of  (-1).  Therefore,  it  follows  that  the  number 
of  equivalence  classes  with  magnitude  |X(a)2,u)2)|  is  at  most  2^P"^^  where  p  is 
the  number  of  non-symmetric  irreducible  factors  in  X(Z2,Z2).  Thus,  the  following 
is  an  immediate  consequence  of  Theorem  4: 

Theorem  5:  Let  xeF(n2,n2)  have  a  z-transform  with  at  most  one 
irreducible  non-symmetric  factor,  i.e. 


X(Z2,Z2)  = 


P(Z2,Z2)  JX 
k=l 
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(52) 


mmmmm 


where  X|^(z^,Z2)  for  k=l,...,p  are  irreducble  symmetric  factors.  If 
ye  F(nj,n2)  with 

I  Y(a)j^,(D2)  1  =  lX(a)j,ui2)  1  (53) 

for  all  0)^  and  1^2*  then  y(nj,n2)  x(nj,n2) . 

As  in  Section  IV,  it  may  be  apparent  that  the  assumption  that  (53)  holds 
for  all  (Dj  and  ui2  in  this  theorem  is  not  necessary  if  x(nj,n2)  and  y(nj,n2)  are 
known  to  be  zero  outside  some  given  domain.  More  specifically,  suppose  that  for 
some  Nj  and  N2,  x(nj,n2)  and  y(n2,n2)  are  known  to  be  zero  outside  the  domain 
^(Nj,N2).  Let  and  be  sets  of  Mj^  distinct  points  as  defined  in  (17)  for 
1(=1,2  and  let  and  .^(Aj.Ag)  be  the  2-D  lattices  generated  by  these 

sets.  Note  that  if 


then 


|X(a)j,a)2) 


=^(^1.^2) 

<^(Aj,A2) 


I  Y(wj,a)2) 


Qy(2l.Z2) 


i?(Ai.A2) 


(54) 

(55) 


where  Qx(Zj»Z2)  and  Qy(Zj,Z2)  are  polynomials,  as  defined  by  (42),  of  degree  at 
most  2(Nj-l)  in  Zj  and  2(N2-1)  in  Z2.  Therefore,  if  Mj>2(Nj-l)  and  M2>2(N2-1) 
then  it  follows  from  Theorem  2  that  Qx^^l’^2^  ” ^y^^l’^2^  ^2* 

(54)  implies  that (53)  holds  for  all  Wj  and  ug.  Consequently, 
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Corollary  5.1:  Let  x,ye  F{nj,n2)  with  support  ^(Nj.Ng).  If  XCzj.Zg) 
has  at  most  one  irreducible  non-symnetric  factor  and 


~  I  X  ((jOj  ) 


(56) 


where  is  a  set  of  distinct  real  numbers  in  the  interval  [0,2ir) 
with  M,^>2(N,^-1)  for  k=l,2,  then  y(nj,n2) '''X(nj,n2) . 


A  special  case  of  this  theorem  results  when  the  points  in  the  sets 
are  equally  spaced  between  0  and  2ir.  In  this  instance. 


lX(m 


1’ 


a>2) 


*  |X(kj,k2) 


Mi,M2 


(57) 


is  the  magnitude  of  the  M^xM2-point  DFT  of  x(nj,n2).  Therefore, (56)  may  be 
replaced  with  the  constraint: 


lY(k^,k2)| 


=  iX(kj,k2) 


M^,M2 


(58) 


4.3  Extensions 


in  Section  4.1,  a  set  of  conditions  are  presented  under  which  a  2-0 
sequence  is  uniquely  defined  to  within  a  scale  factor  by  the  phase  of  its  Fourier 
transform.  A  similiar  set  of  conditions  are  presented  in  Section  4.2  which  allow 
a  2-D  sequence  to  be  uniquely  specified  by  the  magnitude  of  its  Fourier  transform 
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to  within  a  delay,  a  sign,  and  a  180°  rotation.  It  is  of  interest  to  note  that 
these  uniqueness  constraints  are  not  mutually  exclusive.  Specifically,  if 
xeF(nj,n2)  has  a  non-symmetric  irreducible  z-transform,  then  xCn^.ng)  satisfies 
the  constraints  of  Theorems  3  and  S.  Therefore,  the  following  result  is  a  direct 
consequence  of  these  theorems: 

Theorem  6:  If  xeF(ni,n2)  and  has  a  non  symmetric  irreducible 
z-transform,  then  x(nj,n2)  is  uniquely  specified  (in  the  sense  of 
Theorems  3  and  5)  by  either  the  phase  or  magnitude  of  its  Fourier 
transform.  If,  in  addition,  x(nj,n2)  is  known  to  have  support 
^(Nj,N2),  then  the  phase  or  magnitude  of  the  MjXM2-point  DFT  of 
x(nj,n2)  is  sufficient  for  this  unique  specification  provided 
Mj>2(Nj-l)  and  M2>2{N2-1). 

Since  a  non-symmetric  irreducible  z-transform  is  a  sufficient  constraint 
for  a  sequence  xeF(npn2)  to  be  uniquely  defined  by  either  its  magnitude  or 
phase,  it  is  of  interest  to  determine  how  restrictive  this  constraint  is.  In 
other  words,  given  an  arbitrary  2-D  sequence,  is  it  likely  that  its  z-transform 
is  non-symmetric  and  irreducible?  An  answer  to  this  question  may  be  found  by 
considering  the  set  ^^of  all  polynomials  in  two  variables  which  have  degree 
Nj^l  in  Zj  and  N2>1  in  Z2.  With  N  =  (Nj+1)(N2-H) ,  note  that  there  exists  a 
one-to-one  correspondence  between  ^and  R^.  Specifically,  each  polynomial 
pe^  may  be  uniquely  represented  as  a  vector  XpeR^-  Conversely,  each  point 
x^eR*'*  corresponds  to  a  polynomial  pe  .  Now,  consider  the  subset  BSof  ^ 
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which  contains  all  those  polynomials  which  are  reducible.  With  the  correspondence 
noted  above,  represents  a  subset  of  R^.  Since  it  may  be  shown  [9]  that 
^  is  a  set  of  measure  zero  in  "almost  every"  polynomial  in  more  than  one 
variable  is  irreducible.  In  a  probabilistic  setting,  this  result  states  that  a 
polynomial  pe^  is  irreducible  with  probability  one.  It  may  also  be  shown,  in 
a  style  similiar  to  that  in  [9]  ,  that  the  set  of  symmetric  polynomials  represents 
a  set  of  measure  zero  in  Thus,  it  follows  that  "almost  every"  polynomial  in 
two  or  more  variables  satisfies  the  constraints  of  Theorem  6.  Consequently, 

"almost  every"  2-D  sequence  with  finite  support  is  uniquely  defined  (in  the  sense 
of  Theorems  3  and  b)  by  either  the  phase  or  magnitude  of  its  Fourier  transform. 

Although  the  results  which  have  been  presented  thus  far  have  been  confined 
to  sequences  with  finite  support,  an  extension  is  easily  made  to  those  sequences 
whose  convolutional  inverses  have  finite  support.  Specifically,  let  x^(nj,n2) 
denote  the  convolutional  inverse  of  a  2-D  sequence  xCn^.ng),  i.e. 

x(nj,n2)  *  x^.(nj,n2)  =  6(nj,n2)  (59) 

where  6(nj,n2)  is  the  2-D  unit  sample  function.  Now  suppose  that  x(nj,n2)  is 
a  stable  sequence  which  has  a  z-transform  of  the  form 

1 

X(z.  ,ZpI  ®  (60) 

^  ^  P(Zl,Z2) 

where  p(Zj,Z2)  is  a  polynomial  in  ^(z^,Z2).  In  this  case,  the  convolutional 
Inverse  of  x(nj,n2)  has  a  z-transform  given  by 
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X^(Zj,Z2)  =  P(2i*22^ 

so  that  X.  EPlnj.ng).  In  addition,  the  phase  or  magnitude  of  x.(nj,n2)  is 
uniquely  defined  by  the  phase  or  magnitude  of  xin^.ng)  respectively: 


(61) 


lX^(wj.W2)l  =  lX(wj  ,0)2)1 

tlx  (a)j,o)2)  =  -  <(>jj(wj.o)2) 


-1 


(62a) 

(62b) 


Therefore,  if  x(n^,n2)  is  a  stable  sequence  with  a  z-transform  given  by  (60), 
then  x(n.^,n2)  is  uniquely  defined  by  its  phase  or  magnitude  if  the  polynomial 
P(z^,Z2)  satisfies  the  appropriate  constraints. 
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V.  SUMMARY 


In  this  paper,  conditions  have  been  developed  under  which  a  multi -dimensional 
sequence  is  uniquely  specified  by  the  phase  or  magnitude  of  its  Fourier  transform. 
These  conditions  were  shown  to  include  almost  all  sequences  in  two  or  more  unknowns 
which  have  finite  support.  Although  it  was  initially  assumed  that  either  the  phase 
or  magnitude  was  known  for  all  frequencies,  the  uniqueness  constraints  were  then 
extended  to  include  the  case  in  which  the  phase  or  magnitude  was  specified  over 
a  particular  finite  set  of  points. 

The  question  of  uniqueness,  however,  is  only  one  aspect  of  the  more  general 
problem  of  multi -dimensional  signal  reconstruction  from  phase  or  magnitude. 

Although  not  considered  in  this  paper,  another  important  question  concerns  the 
development  of  practical  algorithms  for  reconstructing  a  multi -dimensional 
sequence  from  either  its  phase  or  magnitude.  A  number  of  algorithms  have  been 
developed  for  multi -dimensional  signal  reconstruction  from  phase  which  include 
the  extension  of  1-D  phase-only  reconstruction  algorithms  [6]  to  higher  dimensional 
sequences.  Reconstruction  from  magnitude,  on  the  other  hand,  still  remains  an 
area  of  active  research.  Although  a  few  algorithms  have  been  proposed  (some 
references  may  be  found  in  [?]  ),  none  appear  to  be  suitable  for  the  reconstruction 
of  arbitrary  multi -dimensional  sequences.  Even  beyond  reconstruction  algorithms, 
another  important  issue  which  remains  to  be  investigated  is  the  sensitivity  of 
a  multi -dimensional  sequence  to  measurement  errors  in  its  phase  or  magnitude. 

In  particular,  it  is  of  interest  to  determine  the  effect  on  the  reconstruction 
of  a  sequence  to  noisy  phase  or  magnitude  information. 
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